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In the Israel-Stewart's theory of dissipative hydrodynamics, we have analysed the STAR data on 
(j) meson production in Au+Au collisions at -^5=200 GeV. From a simultaneous fit to <j!> mesons 
multiplicity, mean pr and integrated 112, we obtain a phenomenological estimate of QGP viscosity, 
rj/s = 0.07±0.03±0.14, the first error is due to the experimental uncertainty in STAR measurements, 
the second refiects the uncertainties in initial and final conditions of the fluid. 
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In recent years, there is considerable interest in viscos- 
ity of strongly interacting Quark-Gluon Plasma. String 
theory based models (ADS/CFT) give a lower bound on 
viscosity of any matter 77/s > l/47r 1]. In a perturbative 
QCD, Arnold et al [31 estimated rj/s ~ 1. In a SU(3) 
gauge theory, Meyer [3|] gave the upper bound 77/s <1.0, 
and his best estimate is 77/5=0. 134(33) at T = 1.165Tc. 
At RHIC region, Nakamura and Sakai [3| estimated the 
viscosity of a hot gluon gas as 77/3=0. 1-0.4. Attempts 
have been made to estimate QGP viscosity directly from 
experimental data. Gavin and Abdel-Aziz \d\ proposed 
to measure viscosity from transverse momentum fluctua- 
tions. From the existing data on Au+Au collisions, they 
estimated that QGP viscosity as 77/5=0. 08-0. 30. Exper- 
imental data on elliptic flow has also been used to esti- 
mate QGP viscosity. Elliptic flow scales with eccentricity. 
Departure form the scaling can be understood as due to 
off-equilibrium effect and utilised to estimate viscosity @ 
as, 77/5=0.11-0.19. Experimental observation that ellip- 
tic flow scales with transverse kinetic energy is also used 
to estimate QGP viscosity, 7;/s ~ 0.09 ± 0.015 0, a value 
close to the ADS/CFT bound. From heavy quark energy 
loss, PHENIX collaboration [8] estimated QGP viscosity 
rj/s Ri 0.1-0.16. Recently, from analysis of RHIC data, 
in a viscous hydrodynamics, upper bound to viscosity is 
given 77/s < 0.5 

In the present paper, from a hyd rody namic analysis of 
the recently measured STAR data [ll| on production, 
we have obtained an accurate estimate of QGP viscosity, 
7;/s=0.07 ± 0.03 ± 0.14, the first error corresponding 
the uncertainty in STAR measurements, the 2nd error 
arising from the uncertain initial conditions of the fluid, 
e.g. initial time, initial fluid velocity, freeze-out temper- 
ature etc. As noted in [l2| , several unique features of <j) 
mesons (e.g. hidden strange particle, both hadronic and 
leptonic decay, not affected by resonance decays, mass 
and width are not modified in a medium [13 etc.) make it 
an ideal probe to investigate medium properties in heavy 
ion collisions. For long, strangeness enhancement is con- 
sidered as a signature of QGP formation JJi]. Compared 
to a hadron gas, in QGP, strangeness is enhanced due 



to abundant gg — s- ss reactions. Early produced ss, if 
survive hadronisation can lead to increased production 
of strange particles compared to pp or pA collisions. Ex- 
perimental data do show strangeness enhancement p^ . 
In STAR measurements [ll| also, compared to a pp col- 
lision, in a Au+Au collision, (f> meson production is en- 
hanced. STAR data [ll| also appear to be consistent 
with a model based on recombination of thermal strange 
quarks [1(|. As it will be shown below, STAR data on 

mesons are also consistent with hydrodynamic model 
and are sensitive enough to give an accurate estimate of 
QGP viscosity. 

Space-time evolution of the QGP fluid is obtained by 
solving, Israel-Stewart's theory of 2nd order dissipative 
hydrodynamics. 



d^.T'^" = 



0, 

_J_(7r^''_ 277V<''7i''>) 
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Eq[T] is the conservation equation for the energy- 
momentum tensor, T^^" = (e + p)u^u^ ~ pgt^'^ + tt'"', 
£, p and u being the energy density, pressure and fluid 
velocity respectively, tt'^'^ is the shear stress tensor (we 
have neglected bulk viscosity and heat conduction). Eq[2] 
is the relaxation equation for the shear stress tensor tt^"^ . 
In Eq[21 D = 71^9^ is the convective time derivative, 
V<'^u'^> = ^{y^'u" + V^Tt'') - ^id.u){g>"' - u^u") is a 
symmetric traceless tensor. 77 is the shear viscosity and 
is the relaxation time. It may be mentioned that in 
a conformally symmetric fluid relaxation equation can 
contain additional terms fv^ . 

Eqs lll2l are closed with an equation of state p = 
p{e). Lattice simulations [3 [l^ indicate that the 
confinement-deconfinement transition is a cross over, 
rather than a 1st or 2nd order phase transition. In Figll] 
a recent lattice simulation 49] for the entropy density 
shown. We complement the lattice simulated EOS 
by a hadronic resonance gas (HRG) EOS comprising all 
the resonances below mass 2.5 GeV. In FiglS] the solid 
line is the entropy density of the '"lattice +HRG" ' EOS. 
The entropy density is obtained as. 
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FIG. 1: Black circles are lattice simulation [l9| for entropy 
density. The black line is the model EOS, obtained by para- 
metric representation to the lattice simulations and a hadronic 
resonance gas at low temperature. 



FIG. 2: Viscous fluid (7;/s=0.08) simulation for temporal 
evolution of momentum anisotropy in b=7 fm Au+Au col- 
lision at RHIC. The solid line is the simulation result from 
VISH2+1 and the dashed line is the simulation result 
from AZHYDRO-KOLKATA. 



s = 0.5[1 - tanh{x)]sHRG + 0.5[1 - tanh{x)]siatuce (3) 

with X = Tc^ige MeV, AT = Q.IT^ = IQ.QMeV. 

Compared to lattice simulation, entropy density drops 
slowly at low temperature. 

Assuming boost-invariance, EqsH] and [5] are solved in 
(r = — z^, x, y, rjs — ^Inf^f) coordinates, with a 
code '"AZHYDRO-KOLKATA"',^ developed at the Cy- 
clotron Centre, Kolkata. Details of the code can be 
found in [21| . To show that AZHYDRO-KOLKATA com- 
putes the evolution correctly, in Fig[TJ we have com- 
pared the temporal evolution of momentum anisotropy 
£p = '^T'=^+Tyy> ^ QGP fluid with a calculation of 
Song and Heinz [l7j . Initial conditions are approxi- 
mately same in both the simulations. Within 10% or less, 
AZHYDRO-KOLKATA simulation reproduces Song and 
Hcinz's [l7j result for temporal evolution of momentum 
anisotropy Ep. 

Solution of partial differential equations fEas lll2p re- 
quires initial conditions, e.g. transverse profile of the 
energy density {£{x,y)), fluid velocity {vx{x,y),Vy{x,y)) 
and shear stress tensor {tt'^'^ (x,y)) at the initial time t;. 
One also need to specify the viscosity (r/) and the relax- 
ation time (r^). A freeze-out prescription is also needed 
to convert the information about fluid energy density 
and velocity to particle spectra and compare with ex- 
periment. 

We assumed that the fluid is thermalised at Ti— 0.6 
fm [23] and the initial fluid velocity is zero, Vx{x,y) = 
Vy{x,y) = 0. Initial energy density is assumed to be 
distributed as [13 

e(b, X, y) = e,[Q.lhNpart{h, x, y) + 0.25iVco;/(b, x, y)], 

(4) 



where b is the impact parameter of the collision. Npart 
and Ncoii are the transverse proflle of the average num- 
ber of participants and average number collisions respec- 
tively. Si is a parameter which does not depend on the 
impact parameter of the collision. As will be discussed 
below, we fix it to reproduce experimental data on (f> 
mesons. Finally, the freeze-out was fixed at Tp=150 
MeV [l^l- The shear stress tensor was initialised with 
boost-invariant value, tt^^ = tt*'*' = 2r]/3Ti, 7r^^=0. For 
the relaxation time, we used the Boltzmann estimate 
Ttt = 3ri/2p. 

In Figia for fluid viscosity, Ty/s^O, 0.08, 0.12 and 0.16, 
hydrodynamic predictions for cf) meson multiplicity (^^5^) 

and mean pT (< Px >)? in b=2.3 fm Au-|-Au collisions 
are shown as a function of central energy density. In a 
hydrodynamic model, as expected, both and < pj. > 
increases with increasing energy density, but increase in 
is steeper than in < p^ >. b=2.3 fm Au+Au colli- 
sions corresponds to 0-5% centrality colfisions. In FiglSl 
the blue and red shaded regions represent the STAR mea- 
surements (statistical and systematic error included) on 
(j) mesons multiplicity and mean px in 0-5% Au-|-Au col- 
lisions. One note that irrespective of fluid viscosity, (j) 
meson multiplicity can be fltted in the hydrodynamic 
model by changing the initial energy density, more vis- 
cous fluid requiring less energy density. It is understood. 
In viscous fluid evolution, entropy is generated and fluid 
requires less initial entropy density or energy density. 
However, with increasing 77/ s, mean px also increases, 
even if, (jj meson multiplicity is kept fixed by reducing 
initial energy density. The reason is understood. Initial 
transverse pressure increases with increasing rj/s lead- 
ing to increased < px >■ One note that simultaneous 
fit to STAR data on (p multiplicity and mean px in 0- 



3 



• 


tl/s=0 




■ 


T|/s=0.08 
r|/s=0.12 


0-5% Au+Au 


T 


T|/S=0.16 


Ti=0.6 fm, Tp=150 MeV 



1 .4 




10 


1.3 




g 

>. 




T3 

^ 6 


















4 


1.2 










2 




> 





.] j 














A 


0.10 


1 


1- 
Q. 






V 










0.9 




0.05 


0.8 




0.00 






1.2 


0.7 




> 11 



15 20 25 30 35 40 
EQ (GeV/fm3) 



FIG. 3: 



(color online) Variation of and < p% > with 

central energy density in b=2.3 fm Au+Au collisions in 
AZHYDRO-KOLKATA. The blue and red symbols are the 

predicted and < p% > respectively for ?7/s=0.0, 0.08, 

0.16, 0.2 an3 0.25. The blue and red shaded regions indi- 
cate the STAR measurements (statistical and systematic er- 
ror included) of and < > in 0-5% centrality Au-|-Au 
collisions. 
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FIG. 4: STAR data on the centrality dependence of <\> meson 
(a) multiplicity, (b) integrated vi and (c) mean pt are com- 
pared with hydrodynamical simulation of ideal and viscous 
fluid. 



5% centrality Au-|-Au collisions is obtained only when 
r]l s < 0.12. For higher viscosity, while it is possible to fit 
(j) meson multiplicity, mean px is over predicted. 

TABLE I: The fitted values of the initial central energy den- 
sity (si) and temperature (Ti) of the fluid in b=0 Au-f An col- 
lisions, for different values of viscosity to entropy ratio (rj/s). 

In the last four rows, x^/N for STAR data on < pf, >,V2 
and the combined data sets (^^j— + < Pt >+W2) are shown. 
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In FigHl in three panels, STAR data 11^ on the cen- 
trality dependence of (p meson (a) multiplicity (^^), (b) 

integrated V2 and (c) meanpT (< Pj- >) are shown. Ideal 
or viscous fluid, initialised to fit (jj meson multiplicity in 
0-5% collisions, reproduce (j) meson multiplicity in all the 
centrality ranges of collisions. STAR collaboration mea- 
sured integrated V2 only in 0-5%, 10-40% and 40-80% 
centrality collisions. In mid-central collisions, V2 is re- 
duced by ~ 10%. In very central collisions, elliptic flow 



is very small and ideal and viscous evolution produce 
similar flow. In peripheral collisions, both ideal and vis- 
cous evolutions overestimate the flow. We also observe 
that the STAR data on < > are not explained unless 
r]/s <0.12, as indicated in FiglH 

In tablelH we have tabulated the initial central en- 
ergy density (sj) and temperature (T^) required to fit the 
STAR data on (f) multiplicity. The error in Si is due to the 
statistical+systematic error in STAR measurement. The 
initial energy density of the fluid can be obtained only 
within ~ 10-15% accuracy. In tableHl we have listed the 
~X^ /N for the data sets analysed. Variation of /N of 
the combined data sets with rj/s, is shown in FigjSl In 
FigEl the solid line is a parabolic flt to the combined 

/N , from which we estimate that the best fit to the 
STAR data on the centrality dependence of (j) mesons 
multiplicity, integrated V2 and mean pT is obtained for 
viscosity over entropy ratio as rj/s = 0.07 ± 0.03. We 
have not shown here, but hydrodynamics with rj/s « 0.07 
consistently explain transverse momentum spectra of 
mesons, along with other particles, e.g. pions, kaons. In- 
terestingly, proton data is underestimated in the model 
by a factor of 2. Elliptic fiow data are also reasonably 
well explained. In a later publication detail results will 
be presented. 

The estimate rj/ s — 0.07 ± 0.03 is obtained with fixed 
values of initial time, ri=0.6 fm, freeze-out temperature 
Tf=150 MeV, and initial transverse velocity w, =0. The 
hard scattering contribution to initial energy density was 
also fixed at 25%. The estimate depends on the assumed 
initial and final conditions of the fiuid. In Fig[Bl we have 
studied the dependence of and < > on: (a) ini- 
tial time, (b) freeze-out temperature, (c) hard scatter- 
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ing contribution to initial energy density and (d) initial 
transverse velocity. The fluid was initialised with cen- 
tral energy density £^=29.1 GeV / fm? , corresponding to 
viscosity to entropy ratio 77/5=0. 08. Other conditions 
remaining the same, hydrodynamic predictions for ^^^^ 

and <pi^ > increases respectively by ^ 50% and ^ 10% 
as Ti increases from 0.6 fm to 1.0 fm. If r,; is reduced 
from 0.6 to 0.2 fm, ^ decreases by - 50%, < > 
remains essentially unchanged. For Ti=l fm, STAR data 
on (j) multiplicity can be reproduced if initial energy den- 
sity of the fluid is reduced by ~ 40% (we are assuming 
that dependence of initial energy density is similar 
at ri=0.6 and 1 fm). However, ~ 40% reduction in en- 
ergy density will also reduce < > by ~ 7% (< > 
dependence on initial energy density is weaker than that 
of ^) and STAR data on < > will be over predicted 
by ~ 3%, requiring 55% less viscosity. Arguing sim- 
ilarly, for Ti=0.2 fm, to fit (jj multiplicity, initial energy 
density is to be increased by ~ 40%, which will also in- 
crease < > by ~ 7%, requiring ^ 90% reduction in 
77/s. We estimate systematic uncertainty in 77/s due to 
uncertain initial time (ri = 0.6 ± 0.4 fm) as ^ 90%. 

One of the major source of uncertainty in hydrody- 
namic simulations is the freeze-out temperature {Tp). If 
Tp is increased from 150 to 160 MeV, hydrodynamic pre- 
diction overestimate the STAR data for by 15%, 

< Pj. > remain essentially unchanged (see FigHh). Ini- 
tial energy density of the fluid can be reduced by ^10% 
to fit (j) multiplicity, which will simultaneously reduce 

< > by ~ 2 %. ~ 2% reduction in < > can 
be compensated by increasing rj/s by ~ 7%. If Tp de- 
creases from 150 to 140 MeV, hydrodynamic prediction 
for underestimate the STAR data by ~ 15% and 
overestimate < pj. > by ~ 10%. Initial energy density 
can be increased by ^ 10% to fit ^3^, which will increase 

< > by ~ 2%. ~ 12% increase in < > can not be 
compensated by decreasing -q/s and Tf=140 MeV will 



FIG. 6: (color online) The ratio of hydrodynamic predictions 
for ^^j^ (the black circles) and < > (the red squares) , to 
STAR measurements in 0-5% Au-|-Au collisions as a function 
of (a) initial time (r^), (b) freeze-out temperature (Tp), (c) 
hard scattering contribution to initial energy density, and (d) 
initial transverse velocity {vr = tanh{ar)), are shown. The 
central energy density of the fluid is £^=29.1 GeV/ fm? and 
viscosity to entropy ratio is ?7/s=0.08 . 



be inconsistent with STAR measurements. Indeed, even 
in ideal hydrodynamics, rF=140 MeV overpredict STAR 
data on < >. We estimate the uncertainty in rj/ s due 
to uncertain freeze-out temperature {Tp — 150±10 McV) 
as - 100%. 

In Figint, dependence of and < > on the hard 
scattering contribution to initial energy density is stud- 
ied. If hard scattering contribution is increased from 25% 
to 95%, hydrodynamic predictions for decreases by 

'-^ 20%, but < pi^ > remain essentially unchanged. As 
argued previously, to fit the STAR data on (p multiplic- 
ity, initial energy density is to be increased by ~ 15%, 
which will lead to increase < > by ~ 3%. ^--^3% in- 
crease in pp can be compensated by reducing rj/s by ~ 
55%. If hard scattering contribution decreases from 25% 
to 5%, predicted and < p^ > change marginally 
and estimate of rj/s will not be affected. We estimate 
the uncertainty in rj/s due to uncertain hard scattering 
contribution (0-100%) as --55%. 

Hydrodynamic predictions for of and < pl^ > as a 
function of the initial transverse velocity Vr ~ tanh(ar), 
for a—0, and 0.06 is shown in FiglBJl. Initial transverse 
velocity mainly increase high pp yield and as a increases 
from to 0.06, remain approximately constant, but 

< p^ > increases and the STAR measurement for < pp > 
is overpredicted by ~ 20%. Indeed, non-zero initial trans- 
verse velocity will increase < p^ > requiring lowering of 
77/s. We then estimate uncertainty in 77/s due to un- 
certain initial velocity as ~ 100%. Finite accuracy in 
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the computer code AZHYDRO-KOLKATA, also add to 
the uncertainty in -q/s. In AZHYDRO-KOLKATA, hy- 
drodynamic evolution is computed within ~5% accuracy, 
leading to the ~7% uncertainty in 77/s. Adding all the un- 
certainties is quadrature, systematic uncertainty in 77/s 
is ~ 175%. From the STAR data on (f) meson we then 
estimate QGP viscosity as 77/5=0.0? ± 0.03 ± 0.14, the 
first uncertainty is due to statistical and systematic un- 
certainty in STAR measurements, the second one is due 
to uncertain initial time (Ti=0. 2-1.0 fm), freeze-out tem- 
perature (Tp=140-160 MeV), percentage of hard scatter- 
ing contribution (f=0-95%) and initial transverse velocity 
(a=0-0.6). 

It may be noted that the sources of uncertainties con- 
sidered here is not exhaustive. For example, the un- 
certainty in the initial energy density may not be rep- 
resented in entirety by the Glauber model, by varying 
only the hard scattering fraction. Color Glass condensate 
initial conditions, with larger initial eccentricity, may 
increase the range of 77/s. The uncertainty in freeze- 
out procedure is also not entirely represented by varying 
freeze-out temperature only. Proper treatment of chem- 
ical freeze-out before the kinetic freeze-out, inclusion of 
resonances may also alter the range of uncertainty in rj/ s. 



Also, as mentioned earlier, the relaxation equation Eql2l 
may contain additional terms. While their contribution is 
expected to be smaller than the terms included, estimate 
of viscosity may change if a more complete relaxation 
equation is used. Uncertainty in the initial shear stress 
tensor is also not considered here. We have used boost- 
invariant value as the initial shear stress tensor. While 
over the time scale r^r, Israel-Stewart stress relax to the 
first order value, i.e. to the boost invariant value, it may 
as well be different at the initial time. Range of uncer- 
tainty in "q/s will also increase if uncertainty in initial 
shear stress tensor is taken into consideration. 

The estimate is obtained from experimental data, 
which include the effect of bulk viscosity, if there is any. 
We have neglected bulk viscosity. Neglect of bulk viscos- 
ity, will artificially increase the effect of (shear) viscos- 
ity. In general, bulk viscosity is an order of magnitude 
smaller than shear viscosity. But in QGP, it is possible 
that near the cross-over temperature, bulk viscosity is 
large [H, [23| . Effect of bulk viscosity on particle spectra 
and elliptic flow is studied in [25| . It appears that even 
if small, bulk viscosity can have visible effect on particle 
spectra and elliptic flow. The present estimate then must 
be considered as an upper bound on QGP viscosity. 
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